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Abstract 

This paper introduces and reviews light forces, atom coohng and atom trapping. The emphasis is on the 
physics of the basic processes. In discussing conservative forces the semi-classical dressed states are used rather 
than the usual quantized field dressed states. 

1 Forces 

The idea that light exerts mechanical forces on matter arose early on in astronomy. In 1619 Kepler suggested that 
it was the pressure of sunlight that made comet tails stream away from the Sun (Minogin and Letokhov 1987). And 
he was right. Nearly 250 years later Maxwell's theory of electromagnetism quantified light pressure. The light force 
on individual atoms was observed by Frisch, of Meitner and Frisch fission fame, in 1933. On Earth lasers provide 
the intensities necessary for exerting useful light forces. In the early 1970s Arthur Ashkin at Bell Labs accelerated, 
levitated and trapped micron sized plastic spheres (Ashkin 1972). This work has developed into optical tweezers 
for manipulating small objects. 

Today atoms are routinely cooled and trapped. Light based mirrors, beamsplitters, and lenses for atoms have 
been demonstrated (Adams et al. 1994a). These are the subject matter of atom optics in which atomic de 
Broglie waves are manipulated as in classical optics. Atomic interferometers constructed with atom optics have 
extraordinary sensitivity (Adams et al. 1994b). Promising applications include atomic clocks, lithography, optical 
tweezers, and atom lasers or bosers (Wiseman and Collett 1995). 

In this section we introduce the absorptive and dispersive limits of light forces. Absorptive forces are dissipative, 
and hence useful for cooling, while purely dispersive forces are conservative and hence useful for maintaining 
coherence. Section || deals with cooling: the first step towards trapping. Cooling has two important milestones: 
the doppler limit and the recoil limit. At each "limit" a new approach to cooling must be adopted. At present there 
appears to be no fundamental lower limit to the temperatures which can be achieved. Section ^ surveys methods 
for trapping ions and atoms. For neutral atoms we survey three basic approaches: dissipative traps, conservative 
traps, and magnetic traps. 



1.1 Classical electrodynamics 

The electric field of an electromagnetic wave sets a charged particle oscillating. The Lorentz interaction with 
the wave's magnetic field then pushes the charge in the direction of wave propagation. Since the wave gives 
momentum to the charge it must itself have momentum. Relativistically, the presence of energy and momentum 
in the energy-momentum 4-vector means that field energy in some frame is momentum in another frame. 

A recent experiment at the University of Queensland demonstrated the transfer of field angular momentum to 
micron sized particles (He et al. 1995). The experiment was interesting because it was electromagnetic orbital 
angular momentum rather than photon spin that was transferred (Allen et al. 1992). 

The field's intensity / (W/m^) is given by the modulus of the Poynting vector P = c^€oE_ x S. Using p — Ejt? 
to convert this to a momentum flux gives the radiation pressure (N/m^): 

R = ca^E X S, |i?| = //c. (1) 

The major macroscopic application of this force may well turn out to be solar sailing (Mallove and Matloff 1989). 
The solar flux at the Earth's orbit is about 1400 W/m^, corresponding to a pressure of 4.7 /iN/m^. Perfect reflection 
doubles this pressure. Aluminium coated mylar films a few microns thick have densities of around 5 g/m^ giving 
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a payload free acceleration of around 2 mm/s^. This is enough for interplanetary travel, taking about a year to 
reach Mars. 

The most useful macroscopic light force today is the dispersive force used in optical tweezers. An object 
with higher refractive index than its surroundings is attracted towards high field regions, while one with lower 
index is repelled. This can be understood by considering the momentum changes of light rays refracted through 
the particles. Fig. ^ In section 1.3 we shall derive the corresponding results for atoms. A field detuned above 
resonance, which induces a refractive index n < 1, repels atoms from high field regions, while a below resonance 
field, which induces a refractive index n > 1, attracts atoms (Loudon 1973). 

Fortunately most biological objects have a greater refractive index than water and so can be stably manipulated 
with the waist of a laser beam. Ashkin et al. (1987) showed that live viruses and bacteria could be manipulated with 
microscope based infra-red optical tweezers. The list of biophysical applications to date includes: manipulation of 
single DNA molecules (Chu 1991) and of human gametes (Tadir et al. 1991), cell micro-surgery (Steubing et al. 
1991), and studies of motor proteins (Ashkin et al. 1990) and mitotic chromosomes (Liang et al. 1991). 



1.2 Dissipative forces 

Quantizing the field with plane wave spatial modes exp(ik ■ r) gives field quanta, called photons, with momentum 
hk. The dissipative, or scattering, force is then particularly simple to understand. Each photon absorption gives 
a momentum kick hk to the absorber. Each subsequent spontaneous emission gives a kick of the same magnitude 
but in a random direction, producing no average force. However the resultant momentum fluctuations are what 
limit cooling, as we shall see in section 
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The dissipative force from a laser can slow a counter-propagating beam of atoms. The changing doppler shift 
as the atoms slow can be overcome by chirping the laser frequency or by Zeeman shifting the atomic transition 
frequency with a spatially varying magnetic field. A slowed beam may be transversely cooled by two-dimensional 
optical molasses. This consists of two lasers counter-propagating transversely to the atomic beam and tuned below 
resonance - so called red detuning. Transverse motion doppler shifts the counter-propagating beam towards reso- 
nance and the co-propagating beam away from resonance. Hence more counter-propagating photons are absorbed 
and the atom slows. This process is called doppler cooling. 

There is a limit to the temperature achievable by doppler cooling - the so called doppler limit 

ksTDop = , Doppler limit, (2) 

where 7 is the transition linewidth or inverse lifetime. Note that here and hereafter we use "temperature" in a 
conventional sense to denote the mean energy without implying thermal equilibrium. For the D line of Na, Toop ~ 
240 /iK. The doppler limit is interesting because it depends only on the linewidth. The molasses cools the atoms 
until their kinetic energy equals half the transition energy uncertainty h'y (Stenholm 1986). This is reasonable since 
the atoms are inelastically scattering photons. Scattered photons will typically differ in energy by the linewidth. 
The difference between the scattered photon energy and the absorbed photon energy appears as the atomic kinetic 



energy. In section 2.1 we derive the doppler limit by balancing the cooling rate with the heating rate due to the 



spontaneous emission recoils. 



1.3 Conservative forces & semiclassical dressed states 

Forces arise when there is a spatial gradient in the expectation value of a system's energy, i.e. Hamiltonian. A 
two- level atom interacting with the electric field E_{z, t) has the Hamiltonian (Allen and Eberly 1987), 



H^h{u;A/2){\e){e\-\g){g\)-d-E{L,t) , (3) 

where \g) and |e) are the ground and excited atomic states respectively, huA is the transition energy, and d its 
electric dipole moment operator. Assuming no spatial dependence of loa the Heisenberg equation of motion for the 
momentum p is 

F=^^^[P^H] = -VH = Vid-E{r,t)) , (4) 
dt in — 

since the time derivative of the momentum is the force, F_. Forces arising from a gradient in the Hamiltonian are 
conservative, there is no energy dissipation as in doppler cooling. The force described by equation (^) is known as 
the gradient force or dipole force. 

The eigenstates of the Hamiltonian, equation (^), are called dressed states (Cohen- Tannoudji et al. 1992). The 
energy expectation value is the sum over the dressed states of the product of their energies and populations. Dipole 
forces arise due to gradients in the dressed state energies, or due to dressed state population gradients, or both. 
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In this paper we only consider the semiclassical case in which E_(r_^t) = 5(r) exp(— + f *(r) exp(ia;t) is a 
classical field. Ignoring the kinetic energy and making the rotating wave approximation, the Hamiltonian becomes 

Hsc = n{ujAl2){\e){e\-\g){g\) 

-id{\e){g\£{r)e-'--' - \g){e\£* {r)e^--') , (5) 

where is the field frequency and d is the dipole matrix element related to the transition linewidth by 7 = 
ALo\d'^ / 'ihc^ . We eliminate the time dependence of the field by working in the interaction picture with i7o = 
fi{ijjL/2){\e){e\ — \g)(()\) and with interaction Hamiltonian 

Hscj = -?i(A/2)(|e)(e|-|g)(.g|) 

-zd(|e)(.9|£(r)-|.g)(e|r(r)) , (6) 

where we have defined the field- atom detuning A = — uj^- The semiclassical dressed states are the eigenstates 
of this Hamiltonian. They and their energies are 

|1> =cos(0)|5)+isin(0)|e) , E^^-^Tl ; (7) 

|2> =sin(0)|.g> -zcos(0)|e) , ^ ; (8) 

cos(26i) = -A/7^, sin(26i) = ^7/7^ . (9) 
We have introduced the Rabi frequency and the generalised Rabi frequency 

Vl = 2d\E\l'h, Rabi frequency, (10) 

7?. = \/A^ + il^, generalised Rabi freq. (11) 

The difference between the generalised Rabi frequency and the detuning is referred to as the light shift. It is the 
shift in the transition frequency due to interaction with the light field. 

Consider a ground state atom slowly moving into a laser beam. The adiabatic theorem (Messiah 1966) implies 
that it remains in the Hamiltonian eigenstate continuously connected to the zero field ground state. Which dressed 
state this is depends on the field detuning 

A > ^ cos26l = -1 cos6' = ^ |2) |g) ; (12) 

A < =^ cos26l = 1 ^ cose* = 1 ^ |1) = |g) . (13) 

Consequently the Hamiltonian expectation value is just the eigenvalue ±KIZ/2. In the limit of detuning A large 
compared to the Rabi frequency Q. 

% t, ( 1 / o \ 2\ 

(14) 

For constant detuning only the second term contributes to the gradient corresponding to an effective potential 

This is the mechanical potential produced by intensity gradients of a far detuned field. It has the same sign as 
A = tJL — ^A- So for A > 0, blue detuning, it rises in high intensity regions from which atoms are therefore 
repelled. Otherwise, for A < 0, red detuning, atoms are attracted to high intensity regions. 

This far detuned potential has a simple interpretation in terms of the bare atomic states. The excited state 
population « s/2 is given by equation (|l^), and the atom is excited by absorbing photons with energy ?iA 
greater than that of the excited atom. The energy excess PJi^ is exactly the potential, equation ([l5|); the energy 
excess appears as mechanical potential energy. 

The red detuned dipole force can be used to construct an atom trap called the far off resonance trap or FORT. 
It is essentially optical tweezers for atoms. A FORT trap has been demonstrated to operate at 65 nm detuning, 
confining Rb atoms at 0.4 mK (Miller et al. 1993). 

Since it is non-dissipative the dipole force is most suitable for constructing atom optics devices which preserve 
quantum coherence. This is a requirement for atom interferometry. Proposed atom waveguides constructed from 
hollow optical fibres utilize the dipole force to repel the atoms from the fibre walls (Marksteiner et al. 1994). 
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1.4 A sense of scale 



Let's get an idea of the scale of the forces, masses, accelerations and so on that are typical of atom optics. For 
concreteness we consider a sodium atom which has a mass of M = 23 au = 3.8 x 10~^^ kg. The yellow D-line 
transition is 3Si/2 ^ 3P3/2 with wavelength 589.0 nm, and lifetime 16.4 ns (Adams et al. 1994a). A saturated 
sodium atom has 50% probability to be excited, so on average it absorbs, and emits, one photon every 33 ns or 
3 X 10^ photons/s. The impulse from the absorbed photon momenta is / = (3 x 10^)?ifc = 3.4 x 10"^" (kg m/s)/s. 
The corresponding acceleration of the atom is I/M ft! 10^ m/s^ ~ 10^ g. At this acceleration a 500 m/s atom is 
brought to rest in half a millisecond over about 13 cm. 

The doppler limit temperature, equation (0), of 240 fiK corresponds to a thermal velocity 



«d = \/|^ = 30cm/s . (16) 



By comparison at a temperature of 300 K the average thermal speed in a Maxwellian distribution is ^/SkBTjirM « 
530 m/s. 

The recoil speed from a photon emission or absorption is 

Wr = ^ = 3cm/s . (17) 
M 

A thermal speed of 3 cm/s corresponds to a temperature of about 1 /iK. This recoil limit speed will turn out to 
be the lower limit to speeds attainable by cooling mechanisms for which the cold atoms interact with light. Since 
p = h/X ioi both matter and light the corresponding atomic de Broglie wavelength is, by momentum conservation, 
just that of the photon XdB = Xught = 589 nm. The corresponding kinetic energy Mvl/2 = 1.7 x 10~^^ J. This is 
nearly 10 orders of magnitude smaller than the energy of a photon hu) = 3.4 x 10""'^^ J. It corresponds to a detuning 
of about A = 1.7 X lO^^^/Ti = 1.6 x 10^^ rad/s. 



2 Cooling 

A great variety of cooling methods have been demonstrated (Metcalf and van der Straten 1994). They can be 
classified with respect to two limits representing the lowest temperatures attainable with the particular method: 
the doppler limit and the recoil limit: 

Tdop = ^r-, Doppler Hmit, TRecoU = 77777—, Recoil limit. (18) 
Zkb ZMkb 

For Na these are respectively 240 /iK and 1.2 /iK. Cooling using doppler tuning of an atomic transition is doppler 
limited due to heating by the scattered photon recoils. Sub-doppler cooling utilizes the multi-level nature of atoms 
to approach the recoil limit. It is limited by the fact that the coldest atoms have scattered a last photon and 
hence have at least the photon recoil energy. Sub-recoil cooling requires passive velocity groups for the atoms to 
accumulate in without scattering photons. 

2.1 Doppler cooling 



Doppler cooling was introduced in section 1.2. In this section the doppler cooling limit is derived by balancing 
the viscous cooling with the recoil heating. We assume a two-level atom with velocity v and transition frequency 
oja interacting with ID optical molasses at frequency ujl- The optical molasses consists of two counterpropagating 
laser beams: beam 1 with wavevector k and beam 2 with wavevector —k. The field-atom detuning is denoted 
by A = lol — i^A- From the semiclassical Bloch equations (Allen and Eberly 1987) the steady state excited state 
population is 

1 s 

"~2TT^' (A - fc.i;)2 + 72/4 ' ^ ^ 

where s is the saturation parameter, expressed in terms of the Rabi frequency SI, equation ([lO|), for the atom field 
interaction. 

We now assume that the transition is sufficiently far from saturation s <C 1 that the absorption from each beam 
is independent. Then the force on the atom is the sum of the momentum absorption rate from each beam 

F « hk'jPes - hkjPe,2 

« nk-f^{si- S2), Sj<l . (20) 
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For the small velocities at the doppler cooling limit we can expand s in v = \k.v\/\k\ to get 

/ 2Akv \ , , 

^v^-^o(i±^^:p^j , (21) 

where sq is the zero velocity saturation parameter. Substituting these expansions into the force equation ( pO| ) gives 

2Ak , , 

J^^^fc7^o ^,^^,^^ . . (22) 

For red detuning, A < 0, this is a viscous type damping force proportional to the atomic velocity. The associated 
kinetic energy cooling rate is 

d{Ecooi) ,p , ^, 2Afc 2 

^_ = (F«)=;.fc7^o^j^:^(t') , (23) 

where the angle brackets denote an average over the atomic ensemble. We next consider the photon recoil heating 
of the atom. 

Each spontaneous emission kicks the atom in a random direction with momentum magnitude Tik, producing a 
random walk in momentum space. So after TV emissions the mean square momentum magnitude is (p^) — N{hk)^. 
The total momentum diffusion rate is four times this. One factor of two comes from spontaneous emission from the 
two beams. The other factor of two comes from the absorption, which is also a random process (Cohen- Tannoudji 
1992). Using the spontaneous emission rate 7So/2 we then have the following expression for the kinetic energy 
heating rate 

dt dt 2M'2^'''' ^ ' 

where we have used the zero velocity saturation parameter, as our goal is to find the temperature minimum. At 
equilibrium the sum of this heating rate and the cooling rate, equation (p3|), is zero 

d{Ecool} , d{Eheat) „ ^ ^^_/2\,^ ^ (^.^ 



dt dt A2+72/4' ' M 

1 .A^_n7/2A 7A 

^2 ' 2 2A 8 V 7 2Ay 



(26) 



The last bracketed expression has a minimum of —2 at 2A/7 = —1 A = —27, so this detuning gives the 
minimum mean kinetic energy. This defines the doppler limit energy ksToop/'^, so 

Tdop^^ . (27) 



A result explained heuristically at the end of section 1.2 



2.2 Sisyphus cooling 

The doppler cooling discussed in the previous section only works for light intensities below saturation. Above 
saturation stimulated processes become important and cooling only occurs for blue detuning ! This blue cooling 
is stronger than doppler cooling but does not achieve as low temperatures (Aspect et al. 1986). This is because 
although the damping is stronger so are the fluctuations. 

It is an example of a more general mechanism called Sisyphus cooling after the mythological Sisyphus, king of 
Corinth, who was doomed by Zeus to roll a rock uphill forever. In a standing wave the Rabi frequency f2 varies 
with the spatial mode function from zero at the nodes (zero field) to a maximum at the anti- nodes. Hence so 
does the energy splitting of the semiclassical dressed states hTZ, equation (pi]). Cooling occurs because, as we shall 
see, spontaneous emission is most likely to occur from energy maxima into energy minima. The atom then moves 
up the energy hill and is again most likely to emit at the maximum. On average kinetic energy is converted to 
atom-field interaction energy which is lost by spontaneous emission. 

A detailed explanation using dressed atomic states with a quantized field has been given by Dalibard and 



Cohen- Tannoudji (1985). By contrast we use semiclassical dressed states, which were introduced in section 1.3. 
Any dressed state with an excited state component can spontaneously emit into any dressed state with a ground 
state component. In general a dressed state can spontaneously emit into itself as well as into the other dressed 
state. 

We now explain blue Sisyphus cooling with reference to Fig. ^ In a standing wave field the electric field is zero 
at the nodes and is maximum at the anti- nodes. So the generalised Rabi frequency, equation (flj), varies from |A| 
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at the nodes to VA^ + at the anti-nodes. |2) is the higher energy state and for blue detuning, A > 0, it equals 
the ground state \g) at a node. It increases in energy and mixes in the excited state |e) as the field increases. |1) is 
the lower energy state and for blue detuning it equals |e) at the nodes and it decreases in energy and mixes in the 
ground state with increasing field. Consequently at anti- nodes |2) has its maximum energy and |1) its minimum 
energy while at nodes the opposite holds. But |2) is most likely to spontaneously emit at anti-nodes since there it 
has the largest excited state component. It can decay either to itself, in which case only the photon energy hoj^ 
is lost, or to in which case the larger energy Ti{u!j^ + TZ) is lost. (Recall that the dressed states are defined in 
a picture rotating with frequency ul, so Hcol is the energy zero). In contrast the dressed state |1) is most likely 
to spontaneously emit at nodes, where it is exactly the bare excited state, |e). At a node it has no ground state 
component and hence can only decay to |2), losing energy h(ujL — A). Since this is smaller than h{ujL + TZ) more 
energy than Hlol is lost on average per spontaneous emission. This energy cannot be provided by the photons 
absorbed from the field, which have energy Tiujl- It comes from kinetic energy and hence cools the atom. 

The transfer of kinetic energy to atom-field interaction energy is the Sisyphus mechanism. An atom which has 
just spontaneously emitted is most likely to be in the energy minimum for its dressed state. Any motion increases 
the field-atom interaction energy until it is most likely to spontaneously emit again at the top of the energy hill. 
Thus kinetic energy is transformed into interaction energy which is dissipated by spontaneous emission. 

Consider a sequence of spontaneous emissions which start and end in the same dressed state while passing 
through the other dressed state. The most likely starting point is at an energy maximum and the most likely end 
point an energy minimum. Hence the total process of (at least) two spontaneous emissions dissipates the energy 
difference between the maximum and minimum, which is just the light shift. Sequences which increase the energy 
are possible, but less likely. 

The limiting temperature is determined by the heating due to spontaneous emissions. Emission from |2) is most 
likely to leave the atom in |2) with no cooling effect, but recoil heating still occurs. 

2.3 Spatially dependent optical pumping 

To cool below the doppler limit, equation (|27|), requires multi-level atoms. Since atoms are multi-level and since 
theorists like two-level atoms it is perhaps not surprising that sub-doppler cooling was discovered experimentally 
(Lett et al. 1988). A common feature of the various sub-doppler, but super-recoil, cooling schemes is optical 
pumping. 

Optical pumping is the transfer of population between atomic magnetic sublevels. For definiteness we consider 
a total angular momentum J = l/2<->J = 3/2 transition, Fig. |^. Circularly polarized light is either cr+ or 
(7^ with photon angular momentum +h and —h respectively. An atom in the m = —1/2 ground state sub- level 
can be excited into the m = -1-1/2 excited state sub-level by absorbing a (T+ photon, and into the m — —3/2 
excited state sub-level by absorbing a a~ photon. Similarly the m = -1-1/2 ground sub-level is excited into either 
the m = +3/2 or m ~ —1/2 excited sub- levels by absorbing a or photon respectively. Absorption of 
polarized light increases the magnetic quantum number m while polarized light decreases it. The net result 
of cycles of absorption and spontaneous emission in circularly polarized light is optical pumping: the preferential 
population of ground sub-levels with extremal m. In this case either the m — 1/2 or m = — 1/2 sub-level. In the 
following we are only concerned with weak fields so that excited state population is negligible. 

Spatially dependent optical pumping occurs when the the pumped sub-level depends on position. In 3D this 
is unavoidable and is usually achieved with three pairs of counter-propagating laser beams. Because the atoms 
are moving the state of the atom lags behind the local steady state. Cooling is achieved when this lag is used to 
remove energy from the atoms as they are optically pumped towards local equilibrium. Since polarization gradient 
cooling is one of the more effective of such mechanisms, we consider it next. 

2.4 Polarization gradient cooling 

Polarization gradient cooling uses counter-propagating beams with perpendicular linear polarizations to achieve 
spatially dependent optical pumping, Fig. ^ This is known as the lin perp lin configuration. Atoms lose kinetic 
energy by a Sisyphus mechanism in which the hills are due to light shifts. 3D lin ± lin cooling to within a few 
recoil velocities, about 25 /xK for Na, has been achieved (Cohen- Tannoudji and Phillips 1990). 

The polarization varies with a period of A/2 from linear to (t+, to the orthogonal linear, to a~ , and back to 
linear. Fig. ^. The linear polarization is at angle 7r/4 to the input linear polarizations. This can be seen from the 
expression for the electric field in the lin _L lin configuration 



where x and y are unit vectors in the x and y directions. Using fc = 27r/A and some trigonometric identities gives 



E_ — X sin(fcz — ) -|- y sin(fcz -\- ujt), 



(28) 



E_— {y ~ x) cos(27rz/A) sm{u!t) + {x + y) sin(27rz/A) cos(ajt) 



(29) 
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At z = the field is linearly polarized in the y — x direction. At z = A/8, sin(27rz/A) — cos(27rz/A) = l/\/2 and the 
magnitude of the field is constant in time and hence it is circularly polarized. At z = A/4, sin(27rz/A) = 1 and the 
field is linearly polarized in the x + y direction. At z = 3A/8, sin(27r2;/A) = l/\/2 and cos(27rz/A) = — 1/\/2, and 
the field is again circularly polarized, but this time with the opposite rotation sense. This polarization gradient 
produces spatially dependent optical pumping. 

Lin _L lin cooling uses light shifting of the ground sub-levels. Light shifting was discussed for two-level transitions 
in section 1.3. The total shift of the m = ±1/2 sub-levels is the sum of the shifts due to each of the two transitions 
to which they are coupled by the lasers. The proportion of a~ and cr"*" varies spatially. And since each transition 
has a different coupling strength, as determined by the Clebsch-Gordon coefficients, the light shifts of the ground 
sub-levels depend on position. 

These shifts are negative for red detuning and for the m = —1/2 sub- level are maximum in magnitude where 
the field is polarized, and minimum where it is (t+. The opposite is true for the m = +1/2 sub-level. Fig. ^ At 
the minimum light shift positions optical pumping into the other sub-level is strongest. The result is that atoms 
move uphill until they are optically pumped into the other sub-level. 

This pumping takes energy out of the atom provided the atomic state lags the local equilibrium by somewhat 
less than the optical pumping time Tp. This is the average time for optical pumping to transfer population between 
the ground sub-levels and is proportional to the intensity. Hence lin _L lin cooling works best for velocities Videai 
such that the atom travels A/4 in Tp, Videai ~ '^/(4Tp). Atoms which are moving too fast do not respond to the 
local pumping and hence do not lose energy. Atoms which are moving too slowly are pumped to the other sub-level 
before reaching the top of the hill and hence lose less energy. 

Polarization gradient cooling does not achieve the recoil limit. Cooling closer to the recoil limit can be achieved 
after trapping the atoms in an optical lattice. Adiabatically expanding the lattice cools the atoms. Reduction of 
Cs temperature from the 3 /iK achievable with polarization gradient cooling to 0.7 /iK has been demonstrated with 
adiabatic cooling (Kastberg et al. 1995). This is short of the Cs recoil limit of 0.13 /xK. 

2.5 Velocity selective coherent population trapping - VSCPT 

For cooling below the recoil limit the problem is to avoid the recoil from the last photon emitted from determining 
the temperature. There are two successful approaches: momentum diffusion into a passive zero velocity group, 
and evaporation of the hottest atoms, avoiding photons altogether. Velocity selective coherent population trapping 
(VSCPT) and Raman cooling use the first approach. 

In this section we discuss a simple model of ID VSCPT based on a three-level lambda transition. More detailed 
discussions and the extension to 3D can be found in the 1990 Les Houches lectures of Cohen- Tannoudji (1992). 2D 
VSCPT cooling of metastable He to 250 nK, 16 times below the recoil limit, was reported in 1994 (Lawall 1994). 

We assume that the ground levels are degenerate and that the two counter-propagating lasers have the same 
frequency, and opposite circular polarizations. Fig. ^. The interaction picture Hamiltonian is 

iTvscPT = |^-^A|3)(3|+d£[e*^--|3)(l|+e-^^-|3)(2|+h.c.] , (30) 

where for simplicity we have assumed that the two transitions have the same dipole moment, d. Note that the 
plane wave spatial dependence of the mode in the direction of propagation z has been made explicit. A set of states 
closed under this Hamiltonian is {|NC,p), |C,p), |3,p)}. The first two states are defined by 

\-HC,p)^l={\i,p-nk)-\2,p + nk)) , (31) 

\C,p) = ^{\l,p-hk) + %p + hk)) . (32) 

The "NC" stands for Non-Coupling and the "C" for Coupling. The interesting thing is the action of the Hamiltonian 
equation (|3^) on |NC,p). Since exp(iA:z) generates a momentum space displacement of the momentum eigenstate 

\p) to \p + hk), 

Hvscpt\^C,p) = -L|i^^|l,p_;,fc)-i^i±^|2,p + M)| 

^ + M!\|NC,p)-^|C,p) . (33) 

For p = it is an eigenstate of the kinetic energy and has no coupling to any other state. This is called a dark 
state because it does not interact with the field. VSCPT works by accumulating atoms in the non-coupling states 
|NC,p ~ 0). The non-coupling states are populated via momentum diffusion due to spontaneous emission from the 
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excited state. The coupling amplitude out of them, into \C,p) states is proportional to p and hence their lifetime 
increases as their momentum p decreases. 

In principle there is no minimum temperature for VSCPT. The range of populated non-coupling states near 
p = decreases with time as they are pumped out via the coupling state. This is because the momentum diffusion 
into a particular momentum group is independent of its momentum, while the |NC,p) <-*■ \C,p) coupling amplitude 
is proportional to p, and hence vanishes as p 0. 

The dark state |NC,p = 0) is a quantum mechanical superposition of the atom moving to the left and to 
the right with the recoil momentum hk. In the 2D case the corresponding dark state is a superposition of four 
atomic momentum states. These states were observed in a Paris experiment (Lawall et al. 1994). To confirm the 
superposition, rather than an incoherent mixture, would require interferometry, which remains to be done. 



2.6 Raman cooling 

Like VSCPT Raman cooling works by isolating the atoms near zero velocity from the cooling fields. However there 
is no dark state in Raman cooling, rather the extremely narrow linewidths possible for Raman transitions enable 
tuning of the Raman pulses so that they do not affect the zero velocity group. 

Because the energy difference between the levels in a Raman transition can be very small they can have narrow 
linewidths. This can be used to select out correspondingly narrow velocity groups, and push them around. Kasevich 
and Chu (1992) used this method to push atoms towards zero velocity and achieved 100 nK for Na in ID, which 
is 10 times below recoil. In 2D and 3D they did not quite get to recoil although they beat polarization gradient 
cooling by a factor of nearly 20 (Davidson et al. 1994). 

Raman cooling uses a sequence of Raman pulses from counter-propagating lasers to push atoms towards zero 
velocity. Each Raman transition changes the atomic velocity by 2hk. The pulses are designed so zero velocity 
atoms are unaffected, and the accuracy with which this can be achieved is one of the limits to the cooling. Since 
many Raman cycles are needed for cooling, the atoms must be repumped to their original state. The zero velocity 
state is populated randomly by atoms which have zero velocity after repumping. 



2.7 Evaporative cooling 

Evaporative cooling has achieved the lowest 3D temperatures to date, 170 nK for Rb (Anderson et al. 1995), and 100 
nK for Li (Bradley et al. 1995). These temperatures should be compared to the recoil limit temperatures, equation 
(p^, of Tfiecoii = 230 nK for Rb and TuecoU = 5.5 /iK for Li. This is what made Bose-Einstein condensation of these 
atoms possible. The cooling is achieved by bleeding off the highest energy atoms from the trap. Thermalisation 
of the remaining atoms by elastic collisions produces cooling. After a subsequent adiabatic expansion Anderson et 
al. (1995) reported a temperature of 20 nK. 



The trap is magnetic, section |3.4| , so that no optical heating occurs. Only atoms in magnetic states such that 
their magnetic dipole is attracted to the trap centre are trapped. If the magnetic state is changed the atoms are 
no longer trapped and they escape. Such spin flips are induced by RF fields tuned to a Zeeman splitting frequency. 
Since this frequency depends on the local magnetic field strength, spin flips can be selectively induced in the outer 
parts of the traps. These are precisely the regions occupied by the high energy atoms. Progressive cooling is 
achieved by ramping down the RF frequency. This moves the resonance region towards the centre of the trap, 
where fields are lowest, and allows colder atoms to evaporate. 



3 Trapping 

The Optical Earnshaw theorem states that atoms for which the dissipative force, section |l.2| , is proportional to the 
intensity cannot be trapped by static configurations of laser beams (Ashkin and Gordon 1983). It is true for the 
same reason as the classical Earnshaw theorem is true, section |3.3| . In a region of space without sources or sinks 
for light there is as much energy fiowing in as out, and the dissipative force is in the same direction as the energy 
flow. 

Fortunately the conditions of the theorem can be violated by adding external fields (Pritchard et al. 1986). In 
the case of the highly successful magneto-optical trap this is an inhomogeneous magnetic field which Zeeman shifts 
the transition so that the force is not just proportional to the local intensity. 



3.1 Optical molasses 



Optical molasses, discussed in section 2.1, does not trap atoms. However because the atomic motion is diffusive 
rather than ballistic the atoms can take seconds to travel centimeters (Chu et al. 1985). The dynamics is that of 
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Einstein/Langevin Brownian motion (Gardiner 1985). Using the damping constant of equation (^), 

"^^^^"" A2+^2/4 ' (34) 
we have the Langevin equation for the ID position a; of a doppler cooled atom 

M^^av + N , (35) 

where we have introduced the zero mean noise term N. An equation for (a;^) may be found after multiplying 
equation (||) by x and using {Mv^/2,) = kBT/2 

& = ?Ml + Cexp{-at/M) , (36) 
at a 

where C is an integration constant. For long times the exponential on the right hand side can be ignored and then 

ix')-i4)^^t . (37) 
a 

For Na the root-mean-square displacement over one second is about a few centimeters for T — 100 fiK. 



3.2 Magneto-optical traps: MOTs 

Magneto-optical traps, or MOTs, are now quite common (Wieman et al. 1995). They are also known as Zeeman- 
optical traps, or ZOTs. They use an inhomogeneous magnetic field to produce a position dependent scattering 
force. The magnetic field Zeeman shifts transitions towards or away from resonance with red-detuned, counter- 
propagating, oppositely circularly polarized laser beams. Consider aJ = O^J=l system. The two transitions 
coupling to circularly polarized light are shifted by positive and negative magnetic fields as shown in Fig. ^. The 
magnetic field generated by anti-Hclmholtz coils is also shown in Fig. ||. For an atom displaced towards the a~ 
polarized beam the magnetic field becomes more positive and it Zeeman shifts the a~ transition into resonance 
and the tr"*" transition away from resonance. Hence the scattering from the cr~ beam dominates and the atom is 
pushed back towards the centre of the trap. The magnetic field and laser polarizations are such that atoms are 
pushed back to the centre after displacement in any direction. 

The MOT lasers cool the atoms too. An atom stopping at its turning point is accelerated back towards the 
centre. However once the doppler shift kv exceeds the linewidth 7 it is no longer accelerated. So the maximum 
rebound speed is about 7/fc, corresponding to a temperature of about M{'-f / k)^ / (Metcalf and van der Straten 
1994). This is much higher than the doppler limit, equation (p7|). 



3.3 Ion traps 

Ion traps have certain advantages over optical traps (Blatt 1992). Most importantly the trapping does not rely 
on interaction with light. The potentially extremely small dissipation in ion traps has allowed the demonstration 
of quantum logic gates (Monroe et al. 1995). They are one of the more promising ways of realizing quantum 
computers (Cirac and Zoller 1995). 

Since the divergence of the electric field in free space is zero any region must have field lines both going in and 
coming out. Equivalently there can be no maxima or minima of the electrostatic potential in free space. This is 
the content of Earnshaw's electrostatic theorem. Consequently an ion cannot be trapped by electrostatic fields 
alone. This problem has been solved in two ways: by adding a magnetic field, the Penning trap, or by oscillating 
the electric field, the Paul trap. A quadrupole potential 

V = Vo{r^~2z^) , (38) 

where r is the radial coordinate, is produced by the ring and endcap electrode configuration shown in Fig. 0. 

The Penning trap achieves radial confinement with an axial magnetic field, in the z direction. This inhibits 
radial motion by converting it to circular cyclotron motion about the magnetic field lines. The overall motion of 
the ion is the combination of a slow magnetron rotation ( w 10 kHz) about the z axis and the fast cyclotron motion 
( « 1000 kHz) about the local magnetic field lines. The electrostatic field provides confinement in the z direction. 
The axial oscillation frequency lies between these two. 

The Paul trap adds RF time dependence to the potential equation ( |38| ) of the form Vq = Vdc + Vac cosut. This 
allows a dynamic stability. Again the ion motion has two components; slow approximately harmonic motion in the 
axial and radial directions, and rapid micromotion oscillation about the local position at the driving frequency ( 
w 10 MHz). 
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3.4 Magnetic traps 



Magnetic traps confine atoms by exerting a force on the atomic magnetic dipole moment. Many different trap 
configurations have been proposed (Bergeman et al. 1987). Since they avoid laser heating, and the Coulomb 
repulsion of ions, magnetic traps have been used to Bose-Einstein condense atoms. 

The force is the gradient of the dipole energy — /i • B, where fi is the dipole moment of the magnetic sub-level of 
the atom and B is the magnetic field. For sufficiently slow atomic motion the the magnetic sub-level adiabatically 
follows the changing direction of the magnetic field as it moves around the trap. For example an atom in the local 
ground state remains in the local ground state. Consequently its magnetic moment changes. Atoms in magnetic 
sub-levels whose Zeeman energy — /x-B decreases with magnetic field strength are forced towards the field minimum. 

Anderson et al. (1995) used a variation on the simple quadrupole trap formed with anti-Helmholtz coils (Midgall 
et al. 1985), called the time-averaged orbiting potential, or TOP, trap. This overcomes a problem with the simple 
quadrupole trap caused by its minimum being zero magnetic field. Close to B = the Zeeman levels are nearly 
degenerate and the atomic state may no longer adiabatically follow the magnetic field. Non-adiabatic transitions 
to non-trapped magnetic sublevels can occur, leading to trap loss. 

This problem is particularly bad in quadrupole traps because the field increases linearly in all directions away 
from the minimum, so the field derivative is not zero at the field zero. The TOP trap adds an RF field which 
rotates the field zero in a circle. The time averaged field then has a non-zero, quadratic minimum where the zero 
previously was. Atoms are still lost from the instantaneous zero, but they are high energy atoms with respect to 
the time-averaged potential and hence contribute to evaporative cooling. 

Bradley et al. (1995) used an arrangement of six permanent magnets to produce a field with a non-zero minimum 
in which they reported Bose-Einstein condensation of Li. 

3.5 Gravitational traps 

Gravity is an important consideration for cold atoms. The energy gained by a Na atom falling through 1 mm is 
about 4 X 10~^^ J, or more than twenty times the recoil energy. Gravity can be used to trap atoms in a trampoline 
configuration. Ten bounces of about 3 mm height have been observed experimentally off a parabolic evanescent 
wave reflector (Aminoff et al. 1993). This corresponds to a trapping time of about 0.1 s. This result can probably 
be substantially improved by enhancing the evanescent field. Confinement times at least 50 times longer have been 
obtained using two sheets of blue detuned light to form a "V" in which atoms bounce (Davidson et al. 1995). 

Other reflector shapes such as pyramids and cones have been analyzed (Soding et al. 1995; Dowling and Gea- 
Banacloche 1995). Recoil limited Sisyphus cooling mechanisms which utilize gravity can be incorporated. For 
example if the sub-level changes during the bounce the average force on the way out can be less than that on the 
way in (Ovchinnikov et al. 1995), leading to an inelastic, cooling bounce. 
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Figure 1: Optical ray diagrams showing the origin of dispersive forces on transparent macroscopic particles. The 
arrows represent light rays with thickness proportional to intensity. The columns of arrows represent the light 
beam. The left and right hand circles respectively represent particles with refractive index higher and lower than 
that of the surroundings. Since the light has momentum and it changes direction by refraction, momentum is 
transferred to the particles. The net momentum transfer is determined by the most intense rays, closest to the 
beam centre. 



Figure 2: Mechanism for blue Sisyphus cooling. The top CTirvc is the standing wave field intensity as a function of 
position. The lower two curves are the energies Ei and E2 of the dressed states |1) and |2). The vertical arrows 
denote spontaneous emissions, which are shown at their most probable locations. The diagonal arrows denote the 
atomic motion, e is a number <^ 1. 



Figure 3: Optical pumping inaJ = l/2<-^ J = 3/2 transition. States are labelled by their magnetic quantum 
number, m. The thick double arrowed lines represent transitions induced by either C7+ or a~ polarized light. The 
thin single arrows represent spontaneous emission. 



Figure 4: Mechanism for polarization gradient lin _L lin cooling. A J = 1/2 ground level, as shown in Fig. 3, 
is assumed. The top atomic energy level diagrams indicate the local direction of optical pumping. The curves 
show the energies of the dressed ground sublevels for red detuned light. The pairs of vertical arrows show optical 
pumping at the positions where it is most probable. The horizontal axis shows how the polarization of the field 
varies over half a wavelength due to the perpendicularly polarized counter-propagating lasers. 



Figure 5: Schematic diagram of the atomic levels, lasers, and transitions used in VSCPT. 



Figure 6: The magneto-optical trap. The left part of the diagram schematically shows Zeeman shifting of a 
J = J = 1 transition associated with positive and negative magnetic fields. The right part schematically 
shows the magnetic field and laser configuration in a MOT. The lasers are labelled by their circular polarization, 
either 0-+ or a~ . The resultant net force on atoms at two locations are shown. 



Figure 7: Ion trap electrodes. Only half the ring electrode is shown. The z direction is vertical and the radial 
direction is horizontal. 
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